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In this paper a mean-field theory for the spin-liquid paramagnetic non-superconducting phase of 
the p- and n-type High-T c cuprates is developed. This theory applied to the effective t — t' — t" — J* 
model with the ab initio calculated parameters and with the three-site correlated hoppings. The 
static spin-spin and kinematic correlation functions beyond Hubbard-I approximation are calculated 
self-consistently. The evolution of the Fermi surface and band dispersion is obtained for the wide 
range of doping concentrations x. For p-type systems the three different types of behavior are 
found and the transitions between these types are accompanied by the changes in the Fermi surface 
topology. Thus a quantum phase transitions take place at x = 0.15 and at x = 0.23. Due to the 
different Fermi surface topology we found for n-type cuprates only one quantum critical concentra- 
tion, x = 0.2. The calculated doping dependence of the nodal Fermi velocity and the effective mass 
are in good agreement with the experimental data. 
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I. INTRODUCTION 

Discovered almost 20 years ago, High-T c copper oxides 
still remain a challenge of the modern condensed matter 
physics. It is not only due to unconventional supercon- 
ducting state with a highest superconducting transition 
temperature T c ever observed. Also they reveal evolution 
from an undoped antifcrromagnctic (AFM) insulator to 
an almost conventional, though highly corrclateoi, Fermi 
liquid system at the ovcrdoped side of the phase diagram. 
Between these two regimes the system exhibits strongly 
correlated, or so-called "pseudogap" metallic behavior up 
to an optimal doping concentration x opt ~ 0.16. 

Recent significant improvements of experimental tech- 
niques, especially of the Angle-Resolved Photoemission 
Spectroscopy (ARPES) and the Scanning Tunneling Mi- 
croscopy (STM), revealed many exciting features of this 
doping dependent evolution. First of all, the Fermi 
surface (FS) at low doping concentrations x has been 
measured^. Together with the previous measurements 
on optimally and ovcrdoped samples (see e.g. review^ 
and references therein) , these observations provide a uni- 
fied picture of the doping dependent FS, which changes 
from the "Fermi arcs'— in the underdoped compounds 
to the "large" Fermi surface in the overdoped systems. 
Though this change is smooth, it occurs around the opti- 
mal doping concentration. Also, the observed evolution 
is consistent with the Hall coefficient Rh measurements^. 

The pseudogap behavior observed in ARPES is also 
tracked in the transport measurements. In particular, 
the resistivity curvature mapping over T — x phase dia- 
gram clearly demonstrates crossover between underdoped 
and overdoped regimes^. And the in-plane resistivity p a b 
shows T-linear dependence only around x op t- 

The drastic change of the quasiparticle dynamics 
around the optimal doping was found by the time- 
resolved measurements of the photoinduccd change 



in reflectivity for Bi 2 Sr 2 Cai_ y Dy. y Cu 2 8+< 5 (BSCCO)X. 
Namely, the spectral weight shifts expected for a BCS su- 
perconductor can account for the photoinduced response 
in overdoped, but not underdoped BSCCO. This agrees 
with the observed difference of the low-frequency spectral 
weight transfer in normal and superconducting states on 
under- and over-doped samples 



Meanwhile, the integral characteristics of the system 
demonstrate more smooth behavior upon increase of 
the doping x. The dependence of the chemical poten- 
tial shift A/i on i shows pinning at x < x op t, and 
evolves smoothly at higher doping concentrations^. The 
measured nodal Fermi velocity vp is almost doping- 
independent within experimental error of 20%iii£. The 
experiments involving combination of dc transport and 
infrared spectroscopy revealed an almost constant ef- 
fective electron mass m*/m = 3.8 ± 2 in the under- 
doped and slightly overdoped La^zSr^CuC^ (LSCO) 
and YBa 2 Cu 6 O a (YBCO)^. Low-a; effective mass de- 
pendence contradicts predictions of the Brinkman-Ricc 
metal-insulator transition thcoryii, predicting diver- 
gence of the m* at the point of transition. One of the 
main drawbacks in this theory is that the magnetic cor- 
relations were neglected. Thus the discrepancy with the 
experiment emphasizes the importance of these correla- 
tions in High-T c copper oxides. 

From the theoretical point of view, the description of 
the crossover between almost localized picture and the 
Fermi liquid regime is very difficult. Starting from the 
Fermi liquid approach one may discuss the overdoped 
and, partly, optimally doped region, while for under- 
doped and undoped samples this approach is not appli- 
cable. The strong-coupling Gutzwiller approximation^ 
for the Hubbard model provides a good description for 
the correlated metallic system. This approximation is 
equivalen1ji£ii£ to the mean-field saddle-point solution 
within a slave-boson approach^. At the same time, as 
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shown within 1/d expansion^, with d being the dimen- 
sionality of the lattice, the Gutzwiller approximation is 
equivalent to the case of d — oo. Obviously, for the quasi- 
two-dimensional systems such as High-X^ copper oxides 
this is not a proper limit. The same applies to the Dy- 
namical Mean-Field Theory (DMFT)12i2tt, which is exact 
only for d = oo. In this limit the short-range magnetic 
fluctuations are excluded. It is not a good starting point 
for the system with long-range AFM order at low x and 
short-range AFM correlations in the underdoped region. 

To describe the doping dependent evolution of the low- 
energy excitations we develop a strong-coupling mean- 
field theory for the paramagnetic non-superconductive 
phase of High-T c copper oxides starting from the local 
limit. To go beyond the usual Hubbard-I approxima- 
tion a self-consistcntly calculated static spin-spin and 
kinematic correlation functions are taken into account. 
Within this approximation in the framework of the effec- 
tive t—t'—t" — J* model with ab initio calculated param- 
eters we obtain a doping-dependent evolution of the FS, 
effective mass and nodal Fermi velocity. The analysis of 
the low-energy excitations behavior for p-type cuprates 
yields two quantum phase transitions associated with the 
change of the FS topology. For n-typc cuprates we ob- 
serve only single quantum critical concentration. The key 
aspect in these findings is an adequate description of the 
electron scattering by the short range magnetic fluctua- 
tions accompanied by the three-site correlated hoppings. 

The paper is organized as follows. In Section ITU the 
effective model and the approximations are described. 
The results of the calculations for p-type cuprates are 
presented in Section ITTT1 Also, the critical comparison of 
the t — t'—t" — J* and t — t'—t" — J models is made, 
and the role of short range magnetic order is discussed. 
Section llVl contains results for n-typc cuprates. The last 
Section summarize this study, and the main points are 
discussed. 



II. MODEL AND APPROXIMATIONS 

High-T c cuprates belong to a class of strongly corre- 
lated systems where the standard local density approxi- 
mation (LDA) schemes and the weak-coupling perturba- 
tion theories yield an inappropriate results. To overcome 
this difficulty recently we have developed an LDA+GTB 
method^i. In this method the ab initio LDA calcula- 
tion is used to construct the Wannier functions and to 
obtain the single electron and Coulomb parameters of 
the multiband Hubbard-type model. Within this model 
the electronic structure in the strong correlation regime 
is calculated by the Generalized Tight-Binding (GTB) 
methoc&Sii. The latter combines the exact diagonaliza- 
tion of the model Hamiltonian for a small cluster (unit 
cell) with perturbative treatment of the intercluster hop- 
ping and interactions. For undoped and weakly doped 
LSCO and Nd 2 - K Ce x Cu0 4 (NCCO) this scheme results 
in a charge transfer insulator with a correct value of the 



gap and the dispersion of bands in agreement with the 
experimental ARPES data (see Refi^i for details). 

Then this multiband Hubbard-type Hamiltonian was 
mapped onto an effective low-energy model^i. Parame- 
ters of this effective model were obtained directly from 
the ab initio parameters of the multiband model. The 
low-energy model appears to be the t — t' — t" — J* 
model (t — t' — t" — J model with the three-site corre- 
lated hoppings) for n-type cuprates and the singlet-triplet 
t — t' — t" — J* model for p-type systems. However, for 
x < 0.7 in a phase without a long-range magnetic or- 
der the role of the triplet state and the singlet-triplet 
hybridization is negligible^. Therefore, the triplet state 
could be omitted and in the present paper we will de- 
scribe the low-energy excitations in the single-layer p- 
and n-type cuprates within the t — t' — t" — J* model. 

To write down the model Hamiltonian we use the Hub- 
bard X-operators^: X" <-> Xj' n = \n) (n'\. Here in- 
dex a <-> (n,, n') enumerates quasiparticlc with energy 
u a = £ n (N + 1) — e n '(N), where e„ is the n-th energy 
level of the iV-electron system. The commutation rela- 
tions between X-operators are quite complicated, i.e. two 
operators commute on another operator, not a c-numbcr. 
Nevertheless, depending on the difference of the number 
of fermions in states n and n' it is possible to define 
quasi-Fermi and quasi-Bose type operators in terms of 
obeyed statistics. In this notations the Hamiltonian of 
the t — t' — t" — J* model in the hole representation have 
the form: 

H = 5>o-/i)Xr+ E t f 9 X a i fl X^ 

+ Y, J f9(SfS g -\n f n g )+H 3 . (1) 

Here /i is the chemical potential, Sf is the spin opera- 

tor, s+ = s 7 = sf = \(xp-x**) t 

n f = "^2 Xf' 17 is the number of particles operator, J/ g = 

(7 

2tj g /E c t is the exchange parameter, E ct is the charge- 
transfer gap. In the notations of Ref^i the hopping ma- 
trix elements tf g corresponds to tj g and — tjl for p- and 

n-type cuprates, respectively, and tf g = tj g . Hamilto- 
nian H3 contains the three-site interaction terms: 

TT _ \ ~" tf m t mg / yIT o vSo vOS V crO v Sa V 0a\ 
H 3- ~~ P X ™ X g ~ X f X m X g ) ■ 

(2) 

There is a simple correspondence between X-operators 
and single-electron annihilation operators: a/Xa = 
52 7A<r(oO-Xy i wnere the coefficients J\o-(a) determines 

a 

the partial weight of the quasiparticle a with spin a 
and orbital index A. These coefficients are calculated 
straightforwardly within the GTB scheme. In the con- 
sidered case there is only one quasi-Fermi-type quasipar- 
ticle, a = (0, <t), with 7Ao-(a) = I, and the Hamiltonian 
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in the generalized form in momentum representation is 
given by: 

h = £ (so 

k,<7 k a >P 

+ E E Vpt"' x f ix p< :x r ( 3 ) 

p,(f a, (3, a, a' 

The Fourier transform of the two-time retarded Green 
function in the energy representation, G\(k,E) = 

at ) ) , can be rewritten in terms of the ma- 

k^v I I E 

trix Green function \D(k,E) - 

L J a/3 



xs 



G x {k,E) = Y J 7a. (a)7L (P)D aP (k,E). (4) 



To go beyond the Hubbard-I approximation we have 
to calculate H(k,E). For this purpose we use an equa- 
tions of motion method for the X-operator*^. The exact 
equation of motion for X£ is: 



iXS = [XI H\ =(e -^)XS + L1 

= (s - »)XS + + (£2 - L2(°)). ( 7 ) 



Here £v is the linearized and decoupled in Hubbard-I 
approximation operator L2, 



q(0 ) _ l + a: l + i*E l-i 



2 E r t, 2 



(8) 



The diagram technique for Hubbard X-opcrators 
has been developed^SiSi and the generalized Dyson 
equation^ in the paramagnetic phase ((Xq = (Xq ' ct )) 
reads: 

D(k, E) = \d^(E) P(k, E)t % - P(k, E)V£f (X^) 



Let us define £2 = Ui — L^ ^ and linearize it with re- 
spect to X~: 

k 



\ " rpa/3 v /3 



/-^ k 





T* P X>1 + Lt 



ot{irr) 



(9) 



E(fc,£) P(k,E). 



(5) 



Here, G 1 (£7) is the exact local Green function, 
G^(£) = <5 a/3 / [E - (e n - £„/)], E(fc, £) and P(fc, £) are 
the self-energy and the strength operators, respectively. 
The presence of the strength operator is due to the re- 
distribution of the spectral weight between the Hubbard 
subbands, that is an intrinsic feature of the strongly cor- 
related electron systems. It also should be stressed that 
T,(k,E) in Eq. (JSJ) is the self-energy in X-operators rep- 
resentation and therefore it differs from the self-energy 
entering Dyson equation for the weak coupling perturba- 
tion theory for G\(k, E). 

Within Hubbard-I approximation^ the self-energy 
T,(k,E) is equal to zero and the strength operator 
P(k, E) is replaced by P af3 (k, E) -> P a P = 5 a(3 F a , where 



p / yn,n 

r a(n,n') — \ / 



X 



f /is the occupation factor. 



Taking into account that in the considered paramag- 
netic phase I^X^ 17 ^ — ^X°' ^) = x, with x be- 
ing the doping concentration, after all substitutions and 
treating all fc-independent terms as the chemical poten- 
tial renormalization, the generalized Dyson equation for 
the Hamiltonian (JTJ becomes: 



where T~ 



arc the coefficients of the lin- 



carization. All effects of the finite quasiparticlc lifetime 
are contained in the irreducible part L°^ lrr \ In this pa- 



per we neglect it, 



a(irr) 



0. 



Since the exact equation for the Green function 



is given by E ( ( X5 



iXS 



XL \ ) it is straightforward to find that in our 

k I I E 

nUfi 



approximation T-. corresponds to the self-energy: 



t{k l E) = -t % . 



(10) 



Introducing notations for the static spin-spin correla- 
tion functions 



G 9 - = J2 e' l{f - 9)ir (XfX^) = 2 
f,g 



£e"^ (SZS* ), 

f 

(11) 



and for the kinematic correlation functions 



D(k, E) 



E-(e - ft) 



1 + x 

2 £ 



A> = E e-' i(/ - 9)? (Xf°X /) , 
f,g 



(12) 



x t% 1 — x 



2 Er 



E(fc, A) 



(6) the expression for the quasiparticlc self-energy becomes: 
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Y,{k,E) 



xN ^ 



k- 



l-x 



'^k-q' 



1 



1 + xN 



E 



k — q 



1 — X 



J 4- 





1 + x 2tgt ? 




E c t 


2 -B c t 


k—q 


\ 1 + a: 2 *fe < ~fc-„- 


E c t 


/ 2 


E c t 



(13) 



Here N is the number of vectors in momentum space. 

Until now we have made two major approximations. 
First, we neglected irreducible part of the self-energy, 
iffS lrr \ thus allowing quasiparticlcs to have infinite life- 
time. Second, and this is not so obvious, in Eq. 10 for 
the strength operator P(k, E) we drop out corrections 
beyond Hubbard-I approximation. These dropped cor- 
rections also lead to the finite quasiparticlc lifetime. The 
consequence of these approximations will be discussed 
later. 

Kinematic correlation functions (|12|l are calculated 
straightforwardly with the help of Green function JSJ. 
The spin-spin correlation functions for the t — J model 
with three-site correlated hoppings H 3 were calculated 
in Refi^i. In this paper the equations of motion for the 

spin- spin Green function ^AJ CT \ Xg a ^ were decou- 
pled in the rotationally invariant quantum spin liquid 
phase, similar to Refsi^S. The results of this approach 
for a static magnetic susceptibility for the t — J model arc 
similar to those obtained by other methods^^. Higher- 
order correlation functions appearing due to the H 3 term 
are decoupled in the following way: 

{(x%x?xF I xf)) u - (x**) {{x^xr I x?°)) u 
((x£x?°x°r I xf°)) u - (X[°X^) {{x™ 1 xf)) u 

Thus, higher-order kinematic and spin-spin scattering 
channels are decoupled. 

After the terms proportional to x (2toi/ E c t) being ne- 
glected, the expression for Fourier transform of the spin- 
spin Green function becomes: 



where 




and magnetic excitations spectrum uj$ represented by the 
Eq. (26) of pape»2i. 

The following results were obtained by self-consistent 
calculation of the chemical potential fj., the spin-spin cor- 
relation functions using Green function l|14|) , and the 
kinematic correlation functions (|12[1 using Green func- 
tion |@ with the self-energy 1)13(1. 
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FIG. 1: (color online) The quasiparticle dispersion (on the 
left) and the density of states (DOS, on the right) in the 
paramagnetic phase of p-type cuprate with x — 0.16. The 
position of the chemical potential is denoted by the horizontal 
(green) line. Results within the Hubbard-I approximation are 
shown by dashed (red) and dash-dotted (cyan) curves for the 
model with (t - 1' - 1" - J* model) and without (t-t' - 1" - J 
model) three-site correlated hoppings, respectively. Bold solid 
(blue) curves represent the results for the t— if —if' — J* model 
with the short range magnetic order. 



III. RESULTS FOR p-TYPE CUPRATES 

For LSCO the LDA+GTB calculated parameters are 
(in eV): t = 0.93, t' = -0.12, t" = 0.15, J = 0.295, J' = 
0.003, J" = 0.007. All figures below are in electron rep- 
resentation. 

First of all we would like to stress the important effects 
caused by the three-site correlated hoppings H 3 and the 
renormalizations due to the short range magnetic order. 
Previously, the importance of the three-site correlated 
hoppings in the normal and superconducting phases has 
been demonstrated in Refsi^iiSiiS!. In Fig. ^ we present 
our results for 16% hole doping within different approxi- 
mations. Evidently, introduction of three-site interaction 
terms results in the change of the position of the top of 
the valence band. Therefore, this will become important 
at small x. In AFM phase of the t — J model there is 
a symmetry around (7r/2,7r/2) point. In the paramag- 
netic phase this symmetry is absent. Due to the scatter- 
ing on the short range magnetic fluctuations with AFM 
wave vector Q = (tt,tt) the states near the (tt, tt) point 
are pushed below the Fermi level (see Fignj, thus to- 
tally changing the shape of the FS. In other words, the 
short range magnetic order "tries" to restore the sym- 
metry around (tt/2 1 tt/2) point. In our calculations the 
short range magnetic fluctuations arc taken into account 
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FIG. 2: (color online) Doping dependent evolution of the kine- 
matic (a) and spin-spin (b) correlation functions within the 
t — t' — t" — J* model for p-type cuprates. Index n enumer- 
ates real space vectors connecting neighboring sites: n — 1 
for nearest-neighbors, n — 2 for the next nearest neighbors, 
and so on. In (c) and (d) the doping dependence of the nodal 
Fermi velocity (in units of eV-a with a being a lattice con- 
stant) and the effective mass are shown. 



via the spin-spin correlation functions 1)11(1. 

Our results for the doping dependence of the kinematic 
and spin-spin correlation functions are shown in Fig. [21 
Note, the kinematic correlation functions K n possess a 
very nontrivial doping dependence. For low concentra- 
tions, x < 0.15, due to the strong magnetic correla- 
tions the hoppings to the nearest and to the next-nearest 
neighbors arc suppressed leading to the small values of 
K\ and if 2, while if 3 is not suppressed. Upon increase 
of the doping concentration above x s» 0.15, magnetic 
correlations decrease considerably and nearest-neighbor 
kinematic correlation function K\ increase. Next major 
change sets at x ps 0.23 when the system possesses al- 
most Fermi liquid behavior: K\ becomes largest of all 
K n 's, while the magnetic correlation functions C n and 
the kinematic correlation function K3 arc strongly sup- 
pressed. 

So, we can clearly define two points of the crossover, 
namely x ps 0.15 and x ps 0.23. The system behavior is 
quite different on the different sides of these points, al- 
though there is no phase transition with symmetry break- 
ing occurs. To understand the nature of these crossovers 
we consider the FS evolution with doping concentration, 
presented in FigEI At low x the FS has the form of the 
hole pockets centered around (±7r/2, ±7t/2) point. Then 
these pockets enlarge and at x = 0.15 all of them merge 
together, forming the two FS contours. Up to x = 0.23 
the FS topologically equivalent to the two concentric cir- 
cles with the central one shrinking toward (0, 0) point. 
For x = 0.23 the central FS contour shrinks to the single 
point and vanishes, leaving one large hole-type FS. 

Apparently, the topology of the FS changes drastically 
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FIG. 3: (color online) Band structure (on the left), density 
of states (in the middle), and Fermi surface (on the right) 
evolution with doping concentration x within the t—t'—t" — J* 
model for p-type cuprates. 



upon doping. In particular, it happens at x = 0.15 and 
at x = 0.23. For the first time the "electronic transition" 
accompanying the change in the FS topology, or the so- 
called Lifshitz transition, was described in Refi^. Now 
such transitions referred as a quantum phase transitions 
with a co-dimension= 1 (see e.g. paper—). Note, when 
the FS topology changes at quantum critical concentra- 
tions x\ = 0.15 and at X2 = 0.23 the density of states 
at the Fermi level also exhibit significant modifications. 
This results in the different behavior of the kinematic 
and magnetic correlation functions on the different sides 
of these crossover points. And the changes in the density 
of states at the Fermi level will also result in the signif- 
icant changes of such observable physical quantities as 
the resistivity and the specific heat. 

Also, from the obtained quasiparticle dispersion we cal- 
culated the doping dependence of the nodal Fermi veloc- 
ity vp and the effective mass m*/m (see Fig. |2Jc) and 
(d)). Nodal Fermi velocity does not show steep variations 
with increase of the doping concentration in agreement 
with the ARPES experimcntsii*^. Effective mass m* 
increase with decreasing x and reveals tendency to the 
localization in the vicinity of the metal-insulator tran- 
sition. But this increase is not very large and overall 
m*/m doping dependence agrees quite well with the ex- 
perimentally observed onei£. Note, the non-monotonic 
doping dependence of both these quantities reflects the 
presence of the critical concentrations Xi and Xi ■ 
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FIG. 4: (color online) The same as in Fig|5| but within the 
model for p-type cuprates without three-site correlated hop- 
pings (t — t' — t" — J model). 



FIG. 6: (color online) The same as in Fig|5| but within the 
t — t' — t" — J* model for n-type cuprates. 
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FIG. 5: (color online) The same as in FigOl but within the 
model for p-type cuprates without three-site correlated hop- 
pings (t — t' — t" — J model). 



To analyze the effect of the three-site hopping term H3 
we also calculated the doping dependence of the band 
structure and FS within the t — t' — t" — J model. The 
behavior of the kinematic and the spin-spin correlation 
functions, presented in Fig. 01 is quite different from that 
of the t — t' — t" — J* model. There is only one quantum 
critical point at x » 0.08 and the effective mass becomes 
very large for x approaching zero. For x > 0.1 the evolu- 
tion of the FS and density of states near the Fermi level 
is smooth, without significant changes (see Fig|SJ). Most 
part of the difference to the t — t' — t" — J* model stems 
from the role of H3 in the energy of states near the (it, ir) 
point, thus determining the topology of the FS and the 
physics at low doping concentrations (see Fignj. 
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FIG. 7: (color online) The same as in Fig|^l but within the 
t — t'—t" — J* model for n-type cuprates. 



IV. RESULTS FOR n-TYPE CUPRATES 

Now let us consider n-type cuprates within the t — 
t' - t" - J* model. For NCCO the LDA+GTB calcu- 
lated parameters are (in cV): t = — 0.50, t' = 0.02, t" 
1). 07. J = 0.195, J' = 0.001, J" = 0.004. 

The obtained doping dependence of the kinematic and 
magnetic correlation functions presented in Fig. El There 
is only one crossover point at x « 0.2. Also, in contrast 
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to the p-type results, the most important kinematic cor- 
relation function on the left of this point is K 2 , rather 
than K3. For x > 0.2 the system demonstrates Fermi 
liquid-like behavior with magnitude of kinematic corre- 
lation function decreasing with the distance, and small 
values of the magnetic correlations. 

The role of the short range magnetic order and three- 
site hopping terms in n-type cupratcs is similar to that of 
p-type. In particular, due to the scattering on the mag- 
netic excitations the states near the (it, it) point pushed 
above the Fermi level, and the local symmetry around 
the {tt/2, 7r/2) points is restored, reminding of the short- 
range AFM fluctuations (see Fig. [TJ) . 

Instead of hole pockets around the (±7r/2, ±7r/2) point 
in p-type, here at low x the electron pockets around 
(±7r, 0) and (0, ±7r) points are present. Upon increase 
of the doping concentration these pockets become larger 
and merge together at x — 0.2. For higher concentrations 
the FS appear to be a large hole-like one, shrinking to- 
ward (tt,tt) point. Therefore, no other changes in the FS 
topology other than at x = 0.2 are present. Referring to 
the same arguments as in previous section, we claim that 
in our approach there is only one quantum critical point 
at x n = 0.2 in the n-type cuprates. Th non-monotonic 
change of the effective mass and the nodal Fermi veloc- 
ity is also present at this concentration, as evident from 
Fig.EIc) and (d). 



V. DISCUSSION AND SUMMARY 

To summarize, we have investigated the doping- 
dependent evolution of the low-energy excitations for p- 
and n-type High-T c cuprates in the regime of strong elec- 
tron correlations within the sequentially derived effective 
model with the ab initio parameters. We show that due 
to the changes of the Fermi surface topology with dop- 
ing the system exhibits drastic change of the low-energy 
physics. Namely, for p-type cuprates there exist two crit- 
ical concentrations, x\ « 0.15 and X2 ~ 0.23. Along 
the different sides of these concentrations the behavior 
of the density of states near the Fermi level, of the kine- 
matic and magnetic correlation functions, of the effective 
mass and the nodal Fermi velocity, is drastically differ- 
ent. This let us speak about crossover, or, taking into 
account the accompanying FS topology changes, about 
quantum phase transitions at these quantum critical con- 
centrations. 

For n-type cuprates due to the specific FS topology we 
obtain only one quantum critical concentration, x n ~ 0.2. 

First of all, we would like to comment on the approxi- 
mations made in this work. Since we use the perturbation 
theory with hopping t and exchange J as small parame- 
ters, appropriate for the strongly correlated regime, the 
real part of the corrections to the results obtained will 
be small to the extent of smallness of the higher powers 
of t/E ct and J/E ct . This will result in the small changes 
of the band dispersion and in the fine details of the FS, 



not changing anything qualitatively 

More concerns give the imaginary part of the neglected 
corrections to the strength operator P(k, E) and to the 
self-energy Y>(k,E) through f J °^ lrr \ Application of the 
equation of motion decoupling method to the Hubbard 
model with finite quasiparticlc lifetime^ reveals that the 
results of the mean-ficld-likc approximation is qualita- 
tively correct. Quantitatively, at low doping the imag- 
inary part of the self-energy leads to the hiding of the 
FS portions above the antiferromagnetic Brillouin zone 
((■7T, 0) — (0,7r) line). This results in Fermi arc rather 
than hole pockets at x < x op t (see Fig. 0. Also, we can 
compare our results to the numerical methods, namely, 
to the exact diagonalization studies^. The quasiparticlc 
dispersion of the t — t' — J model in Ref4i can be consid- 
ered as consisting of two bands. For p-type cuprates in- 
tensities of the spectral peaks corresponding to the band 
situated mostly above the Fermi level (in electron rep- 
resentation) are very low. This band is often called a 
"shadow" band and appears due to the scattering on the 
short range AFM fluctuations. Notably, our band dis- 
persion from Figs. |3 and [7| reproduce very well the shape 
of the other, "non-shadow" , band. It is this band where 
the most part of the spectral weight is residing, thus de- 
termining most of the low-energy properties, except for 
such subtle effects as a so-called "kink" in dispersion^. 

Also, all renormalizations not included in considera- 
tion will change the values of the critical concentrations 
X\, x%, and x n . Comparing with the results of a more 
rigorous theory in papcs^S, we expect these values to de- 
crease. 

Thus we conclude that our theory captures the most 
important part of the low-energy physics within the con- 
sidered (and justified for cuprates) model. This claim 
is supported by the qualitative agreement with the crit- 
ical concentrations of crossover observed in the trans- 
port experiments^ and in the optical experimenliLSiS, 
and even quantitative agreement of the doping depen- 
dence of the nodal Fermi velocity and of the effective 
mass ii i i2 i i3_ Although we use a simple mean-field the- 
ory, though strong-coupled, the agreement with the ex- 
periments is not surprising since we included all necessary 
for High-T c copper oxides ingredients. Namely, the short 
range magnetic order and three-site correlated hoppings. 
Former is the intrinsic property of the cuprates exhibit- 
ing long range AFM order at low x, while latter results 
from the sequential derivation of the low-energy effective 
model. 
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